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^ ; ABSTRACT 

^ ■ We present the yield and spectra of stable secondary particles (7, e^, z/g, z/g, 

\^ • u^, and Ufj,) of p—p interaction in parameterized formulae to facilitate calculations 

involving them in astronomical envir onments. The formulae are derived from the 



representing A(1232) and the other representing multiple resonances around 1600 



■ up-to-date p—p interaction model by lKamae et al.l (120051 ) . which incorporates the 

logarithmically rising inelastic cross section, the diffraction dissociation process, 

00 ■ and the Feynman scaling violation. To improve fidelity to experimental data 

l/-^ ■ in lower energies, two baryon resonance contributions have been added: one 

O 

O ■ MeV/c^. The parametrized formulae predict that all secondary particle spectra 

^ ■ be harder by about 0.05 in power-law indices than that of the incident proton and 

their inclusive cross-sections be larger than those predicted hj p — p interaction 
2 ■ models based on the Feynman scaling. 

I Subject headings: cosmic rays — galaxies: jets — gamma-rays: theory — ISM: 

general — neutrinos — supernovae: general 

%. 

1. Introduction 

Gamma-ray emission due to neutral pions produced by proton-proton interaction has 
long been predicted from the Galactic ridge, s upernova remnants (SNRs), active galactic nu- 



cleus (AGN) jets, and other astronomical sites (lHayakawalll969l : ISteckeiill97ll : lMurthy fc Wolfendale 
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19861 : ISchonfelderll200ll : ISchlickeiserll2002l : lAharonianll2004l : lAharonian et al.ll2004f ). High en- 
ergy neutrinos produced by p—p interaction in A GN jets will s oon be detected with large-scale 
neutrino detectors that are under construction (jHalzenll2005l ). Spectra of these gamma-rays, 
neutrinos, and other secondaries depend heavily on the incident proton spectrum, which 
is unknown and needs to be derived, in almost all cases, from the observed spectra them- 
selves. Such analyses often involve iterative calculations with many trial proton spectra. 
The parameterized model presented here is aimed to improve accuracy of such calculations. 

Among the secondaries of p — p interaction in astronomical environment, gamma-rays 
have been best studied. The p-p interaction is one of the two dominant gamma-ray emission 
mechanisms in the sub-GeV to multi-TeV range, the other being Compton up-scattering of 
low energy photons by high energy electrons. Gamma-rays in this energy range ha ve been 
detec t ed from pulsars, the Galactic Ridge, SNRs, blazars and othe r source catego r ies fISteckei 



1971 



20M). 



Murthv fc Wolfendald 1 1989 : lOnd Il998l : ISchonfeldei] I2OOII : IWeekesI l2003l : lAharonian 



High energy gamma-rays from AGN jets are interpreted as mostly due to the inverse 
Compton up-scattering of low-energy photons by multi-TeV electrons. Observed radio and 
X-ray spectra match those of synchrotron radiation by these electrons. Synchronicity in 



(One 


1998; 


Schonfelder 


2001: 



AharonianI |2004| ). For some AGN jets, the above scenario do es not work well, and p — 



p interaction has been proposed as an alte rnative mechanism (IMiicke &: Protheroe 1 12001 
Miicke et allboosl : iBottcher fc Reimeil 12004 ). 



High energy gamma-rays detected by COS-B and EGRET from the Galactic Ridge, on 



tion of protons and nuclei with the interstellar matter (ISM: 


Steckei 19731: Strong et al. 


1978, 


1982. 200c 


): Stephens & 


Badhwai 


1981: 


Maver-Hasselwander et a 


1982|: 


Bloemen et al 


1984; 


Bloemen 


1985: Dermer 


1986a: J 


jteckei 


1990: Hunter et al. 


1997: 


Mori 


1997: Stanev 


2004). 


The measured eamma-rav flux and spectral shape (Hunter et al. 


1997) 


have been viewed as 



the key attestation to this interpretation. It is also known that inver se Compton scattering 



1986; 


Strong; et al. 


2000; 


Schonfelder 


2001) 



In the past two years, several SNRs, including RX J1713-3946 and RX J0852-4622, 
have been imaged in th e TeV band with an angular resolution around 0.1° by H.E.S.S. 
( Aharonian et al. 2004bl fl lAharonian 2005 ). A smooth featureless spectrum, suggestive of 
synchrotron radiation by naulti-T e V electrons, has been detected in the X-r ay band from RX 
J1713-3946 jKovama et al.l Il997l : Isiane et all Il999l : lUchivama eTaP l2003h and RX J0852- 
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4622 (iTsunemi et al.ll2000l : llyudin et al.ll2005l ). The measured TeV gamma-ray fluxes and 
spectra, however, do not a gree well with those pr edicted by the inverse Compton scenario 
(see, e.g., the analysis in Uchiyama et al. ( 20031 )). Several authors have proposed that 
the TeV gamma-rays are possibly due to interaction of accelerated protons with the ISM 
teerezhko fc Volklboool : ISnomoto et al1l2002l : kharonianlbooi iKatadri et aPboOsh . 



Higher precision data are expected, in the GeV range, from GLAST Large Area Tele- 
scope ( GLAST-LAT 2005l[ll an d, in the TeV range, from the upgraded Air Cherenkov Tele- 
scopes ( lAharonian et al.ll2004l ): they will soon test applicability of the inverse Compton up- 
scattering and the proton interaction with ISM for various astronomical gamma-ray sources. 
In many objects, secondary electrons and positrons may produce fluxes of hard X-rays and 
low-energy gamma-rays detectable with high-sensitivity instruments aboard Integral, Swift, 
Suzaku, and NuSTAR^ The formulae given here will give fluxes and spectra of these sec- 
ondary particles for arbitrary incident proton spectrum. 



This work is an extension of that by iKamae et al.l (120051 ). where up-to-date knowl 
edge of the 7r° yield in the p — p inelastic interaction has been used to predict the Galac 
tic d i ffuse gamrna-ray emission. The authors have found that past calculat io ns (IStecker 



197 ^ . 119731 



199G 



Strong et al]ll978l : Istephens fc BadhwailEgSli IPermeijEoSGal lbl: lMorilll997 



Strong et al.ll2000l ) had left out the diffractive interaction and the Feynman scaling violation 
in the non-diffractive inelastic interaction. Another important finding by them is that most 
previous calculations have assumed an energy-independent p — p inelastic cross-section of 
about 24 mb for Tp ^ 10 GeV, whereas recent experimental data have established a log- 
arithmic increase with the incident proton energy. Updating these shortfalls has changed 
the prediction on the gamma-ray spectrum in the GeV band significantly: the gamma-ray 
power-law index is harder than that of the incident proton; and the GeV— TeV gamma-ray 
fiux is s ignificantly larger than that predi cted on the constant cross-section and Feynman 
scaling (IKamae et al.l 120051 ). The model by lKamae et al.l (120051 ) will hereafter be referred to 
as model A. 

Model A does not model p — p interaction accurately near the pion production thresh- 
old. To improve prediction of gamma-rays, electrons, and positrons produced near the pion 
production threshold, two baryon resonance excitation components have bee added to model 
A: A(1232), representing the A resonance, and res(1600), representing resonances around 



^ GLAST Large Area Telescope, |http://www-glast. stanforcl.edu 

^See INTEGRAL Web site, |http://sci.esa.int/esaMI/Integral[ the NuSTAR Web site, 



|http: / /www. nustar.caltech.edu 



the Suzaku Web site 



http: / / www.isas.jaxa.jp/ e /enterp / missions / astro-e2 



and the Swift Web site, http://swift.gsfc.nasa.gov/docs/swift 
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1600 MeV/c^. We note here that A(1232) is the most prominent and hghtest Baryon res- 
onance excited in p — p interaction. It has a mass of 1.232 GeV/c^ and a width of about 



0.12 GeV/c^, and decays to a nucleon (proton or neutron) and a pion 



TT 



-,o,- 



The other 



resonance, res(1600), is assumed to decay to a nucleon and two pions. Introduction of these 
contributions have necessitated adjustment of the model A at lower energies as described 
below. The readjusted model will be referred to as the "readjusted model A". 

The parameterized model presented here exhibits all features of model A at higher 
energies (proton kinetic energy, Tp > 3 GeV) and reproduces experimental data down to the 
pion production threshold. The inclusive gamma-ray and neutrino cross-section formulae 
can be used to predict their yields and spectra for a wide range of incident proton spectrum. 
Formulae for electrons and positrons predict sub-TeV to multi-TeV secondary electrons and 
positrons supplied hj p — p interaction. We note that space-borne experiments such as 
PAMELAlfl will soon measure the electron and positron spectra in the sub-TeV energy range, 
where the secondaries ofp—p i nteraction may become comparable to the primary components 
(lMiillei][200ll : IstephensI EoOll : IPuVernois et al.lboOlh . 



Due to paucity of experimental data and widely accepted modeling, we have not param- 
eterized inclusive secondary cross-sections for a — p nor p—He nor a— He interactions. We 
note that a-par ticles are known to make up about 7% by number of cosmic-rays observed 
near the Earth fjSchlickeiserl |2002| ) and that He to make up ~ 10% by number of interstellar 
gas. The total non p — p contribution is comparable to that oi p — p contribution. The 
a-particle and He nucleus can be regarded, to a good approximation, as four independent 
protons beyond the resonance region (Tp > 3 GeV): th e error introduced i s expected to be 
less than 10% for high energy light secondary particles (IKamae et al.ll2005l ). 



Inclusion of a-particle as projectile and He nuclei as targets will change the positron- 
electron ratio significantly (about 10 — 15%) as discussed below. Fermi motion of nucleons 
and multiple nucleonic interactions in the nucleus are known to significa ntly affect pion 
production near the threshol d and in the resonance region (Tp < 3 GeV; ICrawford et al. 
1980l : iMartensson et al.ll2000l ): we acknowledge need for separate treatment of p— He, a — p 
and a— He interactions in the future. 



See http://wizard.roina2.infn.it/painela 
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Monte Carlo Event Generation 



The parameterization of the inclusive cross sections for 7, e^, z/, and v has been carried 
out, separately, for non-diffractive, diffractive, and resonance-excitation processes, in three 
steps: First, the secondary particle spectra have been extracted out of events generated for 
mono-energetic protons (0.488 GeV < Tp < 512 TeV) based on the readjusted model A. We 
then fit these spectra with a common parameterized function, separately for non-diffractive, 
diffractive and resonance-excitation processes. In the third step, the parameters determined 
for mono-energetic protons are fitted as functions of proton energy, again separately for the 
three processes. The above procedure has been repeated for all secondary particle types. 

The functional formulae often introduce tails extending beyond the energy-momentum 
conservation limits, which may produce artifacts when wide range spectral energy density 
(£'^(ifiux(7)/(i£^) is plotted. To eliminate such artifacts, we introduce another set of functions 
to impose the kinematic limits. 



Several simulation programs have been used in model A ( 



K!amae et al 



high energy non-diffractive process (Tp > 52.6 GeV), Pythia 6.2 (ISiostrand et al. Il200lh with 



the multi-parton-level scaling violation option (ISjostrand fc Skandsll2004f) 



2005): for the 



for the lower en- 



ergy non-diffractive process, the parameterized model by (IBlattnig et al. J2000l ): for the 
diffractive part, the program by T. Kamae (2004, pe^onal communicatio„ii In the .ead- 
justed model A, two programs to simulate two resonance-excitation components have also 
been added. Modeling of the two resonance components will be explained below. 



3. Non-Diffractive, Diffractive, and Resonance-Excitation Cross-Sections 

Experimental data on p — p cross-sections are archived for a broad range of the inci- 
dent proton energy and various final states. T he total and elastic cross-sections have been 
compiled from those by iHagiwara et al.l (120021 ) , as shown in Figure [H The two thin curves 
running through experimental data points in the figure are our eye-ball fits to the total and 
elastic cross sections. We then define the "empirical" inelastic cross section as the difference 
of the two curves to which the sum of non-diffractive, diffractive, A(1232)-excitation, and 
res(1600)-excitation components are constrained. Typical errors in the empirical inelastic 
cross-section are 20% for Tp < 3 GeV and 10% for Tp > 3 GeV. 



*See http://cepa.fnal.gov/CPD/MCTuningl and http://www.phys.ufl.edu/~rfield/cdf 
^Diffractive process has been included in Pythia after the work began. 
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The four component cross sections of the readjusted model A and their sum are shown 
in Figures [1] and [21 The empirical inelastic cross section is shown by a series of small circles in 
Figure[2J The component cross sections take formulae given in equations ([T]) (non-diffractive), 
dg) (diffractive), Q [A(1232)] and (HD [res(1600)]. These are also shown in Figures □ and H 

We note that there is no clear experimental method separating the four components, 
especially at lower energies (Tp < 20 GeV). This ambiguity does not significantly affect the 
secondary particle fluxes, as long as the sum agrees with the total inelastic cross section and 
the total secondary inclusive cross sections agree with the corresponding experimental data. 

The secondary particle spectra for a mono-energetic proton are normalized to the com- 
ponent cross sections given in equations ([I]), ([2]), ([3]) and (jl]) at the corresponding proton 
energy. W e note that the non-di ffractive component for Tp < 52.6 GeV is based on the 
formula by iBlattnig et al. I (120001 ). which is normalized to their 7r° inclusive cross section 
formula not to the total inelastic cross section. In the readjusted model A, this compo- 
nent cross section is defined by equation ([1]) and the 7r° inclusive cross section formula of 
Blattnig et al. I (120001 ) has been redefined so that the sum of the four components reproduce 
the experimental vr^ inclusive cross section. The positive and negative pion inclusive cross 
sections are also redefined as products of ou r tt" inclusive cr o ss sec tion and the ratio of the 
7T~^'~ and TT^ inclusive cross sections given in iBlattnig et al. I (120001 ). 
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Pp [GeV/c] 



Fig. 1. — Experimental p — p cross sections, as a function of proton momentum, and that 
of readjusted model A: experimental total (squares), experimental elastic (triangles), total 
inelastic (thick solid line), non-diffractive (dashed line), diffractive process (dot-dashed line), 
A (1232) (dotted line), and res(1600) (thin solid line). The total inelastic is the sum of 
the four components. The thin solid and dot-dot-dashed hnes running through the two 
experimental data sets are eye-ball fits to the total and elastic cross sections, respectively. 
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12 5 10 



Pp [GeV/c] 

Fig. 2. — Experimental p — p cross sections, as a function of proton momentum, and that 
of readjusted model A for Tp < 10 GeV. Small circles represent the empirical inelastic cross 
section described in the text. Lines are the same as in Figure [H 
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^NonDiff(^)[mb] 



Pp<l GeV/c, 
0.57(a;/ao)^'^(a2 + asx^ + a^x^ 

+05 exp(-a6(x + 07)^)) 1 < -Pp < 1-3 GeV/c, 
(6o|ai -x\+ bi\ao - x\)/{ai - oq) 1.3 < Pp < 2.4GeV/c, 
02 + 03^^ + a4X^ 

+a5 exp(-a6(a; + 07)^) 2.4 < Pp < 10 GeV/c, 

Co + Cix + C2a;^ Pp > 10 GeV/c, 



(1) 



'^mfr(^)[mb] 



r 

^/(a; - (io)/c^i 
x{d2 + d3\ogiQ{d4{x - 0.25)) 

(i2 + d3 logio{d4{x - 0.25)) 
+(i5X^ - dex^ 
[ co + eix 



Pp < 2.25 GeV/c, 



2.25 < Pp < 3.2 GeV/c, 

3.2 < Pp < 100 GeV/c, 
Pp > 100 GeV/c, 



(2) 



'A1232 



l^a;) [mb] 





/iexp(-/2(Ep-/3)2) 




Ep < 1.4 GeV, 
1.4 <Ep< 1.6 GeV, 
1.6 <Ep< 1.8 GeV, 
1.8 < Ep < 10 GeV, 
Ep > 10 GeV, 



(3) 



Res(1600) 



[x) [mb] 





9oE. 



14 



£;p < 1.6 GeV, 
1.6 <Ep< 1.9 GeV, 
gi exp (-^2(£^p - ^3)') 1.9 <Ep< 2.3 GeV, 
94Ep^ 2.3 < ^p < 20 GeV, 

£;p > 20 GeV, 



(4) 



where x — log^Q(Pp[GeV/c]) and Ep is the proton energy in GeV. 



3.1. Introduction of Resonance-Excitation Processes to Model A 



One or both of the projectile and target protons can be excited to baryon resonances in 
the p— p interaction. Here we use "baryon resonances" to represent both nucleon resonances 
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(iso-spin=l/2) and A resonances (iso-spin=3/2). These excitations enhance the pion pro- 
duction (and hence secondary particle production) near the inelastic threshold. The most 
prominent resonance among them is A(1232), wh ich has a mass o f 1232 MeV/c^ and decays 
predominantly (> 99%) to a nucleon and a pion (IHagiwara et al.ll2002l ). 

Steckeil Jl97oh proposed a cosmic gamma-ray model in which neutral pions are produced 



only through the A(1232) excitation for Tp < 2.2 GeV. The resonance is assumed to move 
only in the direction of the incident proton. At higher energies, another process, the fireball 
process, sets in and produces pions with limited transverse momenta. 



Dermerl (jl986al ) compared predictions of models on -k^ kinetic energy distribution in 



th e proton-proton center-of-mass (CM) system with experiments and noted that the model 
bv ISteckerl (Il970l ) reproduces data better than the scaling model by [Stephens &: Badhwar 
( Il98ll ) for Tp < 3 GeV. He proposed a cosmic gamma-ray production model that covers a 
wider energy range by connecting the two models in the energy range Tp = 3 — 7 GeV. 



Model A by iKamae et al.l (120051 ) has been constructed primarily for the p — p inelastic 
interaction Tp ^ 1 GeV and has left room for improvement for Tp < 3 GeV. The diffrac- 
tion dissociatio n compo r ient o f model A has a resonance-excitation feature similar to that 
implemented in ISteckerl (jl97ll ) for Tp > 3 GeV where either or both protons can be excited 
to nucleon resonances (iso-spin=l/2 and mass around 1600 MeV/c^) along the direction of 
the incident and/or target protons. What has not been implemented in model A is the en- 
hancement by baryon resonances in the inclusive pion production cross sections below Tp < 3 
GeV. 



We note h ere that the the models by IStecked (Il970l ) and by iDermerl (Il986al . see also 
Dermerlll986bl ) used experimental data on the inclusive 7r° yield (and that of c harged pions) 



to gui de their modeling, but not the total inelastic cross section. Model A by iKamae et al. 



(120051 ) ■ on the other hand, has simulated all particles in each event (referred to as the "ex- 
clusive" particle distribution) for all component cross sections. One exception i s simu lation 
of the low-energy non- diffract ive process (Tp < 52.6 GeV) by iBlattnig et al. I (120001 ). The 
inclusive 7r° (or gamma-ray) yield is obtained by collecting vr" (or gamma-rays) in simulated 
exclusive events. When r eadju sting model A by adding the resonance-excitation feature 
similar to that by ISteckerl (Il970l ). overall coherence to model A has been kept. We adjusted 
the A(1232) excitation cross section to reproduce the total inelastic cross section given in 
Figure [2] and fixed the average pion multiplicities for + : : — to those expected by the 
one-pion-exchange hypothesis, 0.73 : 0.27 : 0.0. As higher-mass resonances begin to con- 
tribute, the average pion multiplicity is expected to increase. To reproduce the experimental 
77° inclusive cross section and total inelastic cross section for Tp < 3 GeV, we introduced a 
second resonance, res(1600). This resonance does not correspond to any specific resonance 
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but represents several baryon resonances at around 1600 MeV/c^: its pion multiplicities 
(+ : : — ) are assumed to be 1.0:0.8:0.2. We note here that the resonance components favor 
positive pions significantly over neutral pions while negative pions are strongly suppressed. 

The distribution of pion kinetic energy in the p—p center-of-mass ( CM) sys t em (TV ) has 
been adjusted to reproduce the experimental ones given in Figs. 3-5 of iDermerl (jl986al ). For 
the A(1232) excitation, the probability increases proportionally to up to its maximum, set 
at = 0.28xabs(Tp— 0.4)°-*^^. Here and Tp are measured in GeV/c. The distribution goes 
to zero beyond this maximum value. For res(1600), the probability distribution increases 
proportionally to to reach its peak at = 0.16 x abs(Tp — 0.4)°-''^. It decreases linearly 
until reaching zero at twice the peak of T,r- 

Pion momentum is directed isotropically in the p — p CM system for res(1600) as well as 
for A(1232). No angular correlation has been assumed between the two pions from res(1600): 
this is justified for the astronomical environment where chance of detecting two gamma rays 
from a same interaction is null. Decay kinematics including the polarization effect has been 
implemented to the charged pion decay. This treatment allows the resonances to recoil 
transversely to the direction of the incident p roton while the recoil was constrained along 



the incident proton direction in IStecked fll97' 



To validate the resonance components of the readjusted model A, we have compared the 
model vr" spectrum in the p—p CM system at Tp = 0.65, 0.97, and 2.0 GeV with experimental 
data in Figures [3], HI O Shown in these figures are contributions of the A(1232), res(1600), 
non-diffractive, and diffractive processes. Our model reproduces well the shape of pion 
kinetic energy distribution at Tp = 0.65 GeV but begins to concentrate more towards zero 
kinetic energy than experimental data at Tp = 0.97 and 2.0 GeV. Fidelity to the experirn ental 
data is much improved when cor npared with th e model by [Stephens fc Badhwan (Il98ll) but 
somewhat worse than the one by IStecked (jl970l ) given in Figures 2-6 of iDermerl (jl986al ). The 
difference among the models becomes less noticeable for pion decay products, gamma-rays, 
electrons, positrons, and neutrinos. 

Our inclusive tt" c ross sect i on, su m of all four compoii ents, is compared with experimen- 
tal data assembled by ISteckerl (Il970l ) and IDermerl (Il986al ) in Fig. [HI The readjusted model 
A reproduces experimental data quite well for a wide range of incident proton energy. 
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0.5 




T^o in CM [GeV] 



Fig. 3. — Experimental and simulated 7r° kinetic energy distributions in the p — p C M system 
for Tp = 0.65 GeV. The thick solid line is the experimental data taken from Fig. 3 of iDermer 
The dashed line is the sum of all readjusted model A components: non-diffractive 
(dot-dashed line), A(1232) (double-dot-dashed line), and res(1600) (thin solid line). 
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Fig. 4.- 
for T„ = 



- Experimental and simulated vr*^ kinetic energy distributions in the p — p CM system 
0.97 G eV. The thick solid line represents experimental data (taken from Fig. 4 of 



Dermerl (jl986al )). Other lines are the same as in Fig. [31 
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Fig. 5. — Experimental and simulated vr" kinetic energy distribution in the p — p CM system 
for Tp = 2.0 G eV. the thick solid line represents experimental data (taken from Fig. 5 of 
Dermerl (jl986al )). and the dotted line is for the diffractive process. The others are the same 
as in Fig. [31 
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4. Inclusive Spectra of Simulated Events for Mono-Energetic Proton Beam 

The first step of parameterization is to generate simulated events for mono-energetic 
protons. To simplify this step, events have been generated for discrete proton energies at a 
geometrical series of Tp = 1000.0 x 2(^~^^)/^ GeV where = — 40. Each proton kinetic 
energy (Tp) represents a bin covering between 2~^-'^^Tp and 2^-'^^Tp. The sampling density 
has been increased for Tp < 1 GeV by adding points at Tp = 0.58 GeV and 0.82 GeV. 

Secondary particle spectra are histogrammed from these simulated events in energy bins 
of width AE/E = 5 %. Figures [7^ and [7b show thus-obtained inclusive gamma- ray cross 
sections for the non- diffract ive and diffractive processes, for Tp = 512 TeV, 8 TeV, 125 GeV 
and 11.1 GeV respectively. Those for are given for Tp = 8 TeV and 512 TeV in Figures 
|8}i and [8]d. 
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Pp [GeV/c] 



Fig. 6. — Exp erimenta l and s imul ated ti^ inclusive cross section. Experimental data are those 
assembled by ISteckerl (Il970l ) and iDermerl (Il986al ). Lines are those of readjusted model A: 
total (solid line), non-diffractive (dashed line), diffractive process (dot-dashed line), A(1232) 
(double-dot-dashed line), and res(1600) (dotted line). 
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Fig. 7. — Gamma-ray inclusive cross sections by readjusted model A for four mono-energetic 
protons for (a) the non-diffractive process and (b) the total (top four curves) and the diffrac- 
tive process (bottom four curves). Note that the diffractive process produces two humps at 
Tp = 512 TeV, 8 TeV and 125 GeV. The four proton kinetic energies are Tp = 512 TeV (solid 
line), 8 TeV (dashed line), 125 GeV (dot-dashed line) and 11.1 GeV (dotted line). The bin 
width is Alog(ii^) = 0.05. The large fluctuation in the curves at the highest and lowest ends 
are due to statistics in the simulation. 
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Fig. 8. — Electron and positron inclusive cross sections by readjusted model A for Tp = 8 
and 512 TeV for (a) the non-diffractive and (b) diffractive processes. Electron spectra are 
shown by solid lines and positron spectra by dashed lines. Note that the spectra have two 
humps in the diffractive process. The bin width is Alog(£') = 0.05. The large fluctuation in 
the histograms at the highest and lowest ends are due to statistics in the simulation. 
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We then define, after several iterations of fitting, functional formulae tliat reproduce 
tfie secondary particle spectra for mono-energetic protons, for the non-diffr active, diffractive, 
and resonance-excitation processes separately. For the non-diffractive process, the differential 
inclusive cross section (Acjnd) to produce a secondary particle in a bin of width AEsec/Egec = 
100 % centered at Egec is given as 

A(Tnd(-E'scc) 77 / M7 / \ /cN 

Al0g(E3ec) = 

where i?sec[GeV] is the energy of the secondary particle and x = log]^Q(i?sec[GeV]). 
Fmix) is the formula representing the non-diffractive cross section, given in equation 
below, and -Fnd.uI^;) is the formula to approximately enforce the energy- momentum conser- 
vation limits: 



Fmix) = aoexp(-ai(x - as + a2(x - as)^)^) + 

tti exp(— a5(x — as + a^i^x — a^Y + iii.^ — o-s)^)^); (6) 



i^ND,ki(a;) ^^^p {W^^^,{L„,,,, - x)) + 1) (exp {W^^,^{x - L^,,)) + 1) ' 

where Lmin and Lmax are the lower and upper kinematic limits imposed and PVnd,i and 
W^ND,h are the widths of the kinematic cut-offs; Lmin = —2.6 for all secondary particles, and 
the other parameters are listed in Table [2j 

For the diffractive process we use a similar function 



Alog(Escc) 



where i?sec[GeV] is the energy of the secondary particle and x = log]^Q(i?sec[GeV]); 
Fdm{x) represents the diffractive cross section, given in equation Qj below, and -Fki(x) en- 
forces the energy-momentum conservation: 



i^Difr(a;) = 6o exp(-6i((x - 62)/(l + &3(a; - 62)))^) + 



(9) 
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with Wdiflf = 75 and L^ax = log^Q{Tp[GeV]). 

For the resonance-excitation processes [A(1232) and res(1600)] we use the function, 

(11) 

where i?sec[GeV] is the energy of the secondary particle and x = log]^o(-^sec[GeV]), F^es{x) 
represents the cross section, given in equation (fT2l) below, and -Fki(x), which is the same as 
for the diffraction process, enforces the energy-momentum conservation: 



Fj.es{x) = coexp (-ci((x - C2)/(l + C3{X - C2) + C4{x - 02)^))^). (12) 

To ensure that the parameterized model reproduces the experimental 7r° multiplicity 
after the readjustment in the resonance-excitation region of Tp, we have renormalized the 
non-diffractive contribution by multiplying it with a renormalization factor, r{Tp), given 
below, to the final spectrum. Note that this readjustment does not affected the diffractive 
process. 



r(Tp) ~ 1.01 for Tp > 1.95GeV, (13) 



r(y = logio(Tp)) =3.05exp(-107((y + 3.25)/(l + 8.08(?/ + 3.25)))2) for Tp < 1.95GeV. 

(14) 

For all other secondary particles, r{y = \og^Q{Tp)) is found in Tables miHl 

The simulated gamma-ray inclusive cross sections are superimposed with the parame- 
terized ones in Figures [9^ and [9]d for three mono-energetic protons, for the non-diffractive 
and diffractive processes and for the resonance-excitation processes in Figure [lOl The agree- 
ment is generally good except near the higher and lower kinematical limits where we find a 
difference of about 10-20%. 
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[GeV] 




10-2 1 102 10-1 



E., [GeV] 

Fig. 9. — Simulated (histograms) and parameterized (solid lines) gamma-ray inclusive cross 
sections for (a) non-diffractive and (b) diffractive processes for Tp = 8, 64 and 512 TeV. The 
parameterized cross sections are defined by equations ([6]) and ([9]), and parameters given in 
Table El 
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Fig. 10. — Simulated (histograms) and parameterized (lines) gamma-ray inclusive cross 
sections for the resonance-excitation processes at Tp = 0.69 GeV. A(1232) (double-dot- 
dashed line), res(1600) (thin solid line), and the sum of all model distributions (dashed 
line). The parameterized cross sections are defined by equations (fT2l) for both A(1232) and 
res(1600), and parameters are given in Table [3l 
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5. Representation of Parameters as Functions of Incident Proton Energy 

The parameterization formulae for secondary particles for mono-energetic protons (equa- 
tions [6], [9] and [12]) have nine, eight, and five parameters for each Tp for non- diffract ive, 
diffractive, and resonance-excitation processes respectively. These parameters depend on the 
proton kinetic energy, Tp. The final step of the parameterization is to find simple functions 
representing energy dependence of these parameters. Functions obtained by fitting often 
give values different significantly from those found for mono-energetic protons near the kine- 
matic limits and produce artifacts in the wide range spectral energy density, as described 
previously. Some manual adjustments have been made to control possible artifacts. 

5.1. Parameterized Gamma Ray Spectrum 

The final functional representation of inclusive cross sections for secondary gamma-rays 
is given in equations (Q, ([9]), and (fT2|) with parameters defined as functions of Tp in TeV 
(not GeV) in Table [31 The total inclusive gamma-ray spectrum is the sum of the non- 
diffractive, diffractive, and resonance-excitation contributions. The spectrum produced by 
protons with a continuous spectrum can be calculated by summing over the total gamma- 
ray spectra for mono-energetic protons with appropriate spectral weight. For example, the 
spectra for power- law protons extending to Tp = 512 TeV with index=2 and index=2.7 
have been calculated and compared with the corresponding histograms produced from the 
simulated events in Figure [TTJ The parameterized model reproduces either spectrum within 
10%: it predicts 10-20% more gamma-rays than simulation by the readjusted model A at 
the higher kinematical limit. 
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Fig. 11. — Gamma-ray spectra produced by protons with power-law spectrum cutoff at 
Tp = 512 TeV: (a) index=2 and (b) index=2.7 by the parametrized model [solid lines) defined 
by equations ([6]), ([9]), f|T2|) . and by Table [31 The histograms are those of the simulation by 
the readjusted model A. The dashed straight line corresponds to an index of (a) 1.95 and 
(b) 2.68. 
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5.2. Parameterized e and Neutrino Spectra 

The parameterization lias been extended to otlier secondary particles, e^, e"*", ^e-, 
Ufj^, and z/^, in tlie same way as for gamma-rays. Tlieir functional formulae are represented 
by equations IQ, ([9]), and f|T2|) with the parameters defined in Tables IH E], O O [H and [9l 
respectively. We note that no is produced in A(1232) decay in readjusted model A and, 
hence, no e~ and Ue either. 

We note that the secondary electron and positron spectra from charged pion and muon 
decays have been calculated including the polarization effect of the weak interaction theory. 
The spectra produced by power-law protons of index=2.0 and 2.7 (Tp < 512 TeV) have been 
computed based on these parameterized models in Figure [12] for and Figure [13] for u^. 
We note in Figure [12] that more are produced than e~ throughout their spectra. This 
is largely due to the charge conservation and enhanced by the fact that we have neglected 
a-particles and neutron decays. The number of electrons produced in the p — p interaction 
will match that of positrons if we include electrons coming out of neutron decays. 

For given Tp, electrons from neutron decays have low energy (mostly with E < 10 MeV) 
as shown by the lower histograms in Figure [121 They do not contribute to the high energy 
gamma-ray spectrum. 
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(b) 




E,± [GeV] 



Fig. 12. — Electron and positron spectra produced by protons with power-law spectrum 
cutoff at Tp = 512 TeV: (a) index=2 and (b) index=2.7 by the readjusted model A and by 
the parameterized model (solid lines) defined by equations IQ, ([9]), and ( fT2ll . The param- 
eters are defined in Table H] for electrons and in Table [5] for positrons. The histograms are 
those of the simulation by the readjusted model A. The dashed straight lines correspond to 
asymptotic power-laws: for electrons 1.93 and 2.63 respectively; and for positrons 1.96 and 
2.69 respectively. The histograms in the lower left corner of (a) and (b) show the spectra of 
the electrons from the neutrons produced by the power-law protons. 
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Fig. 13. — Muon neutrino spectra produced by protons with power-law spectrum of (a) 
index=2 and (b) index=2.7 by the parameterized model (solid curves) defined by equations 
([6]), ([9]), and (fT2ll . The parameters are defined in Table [HI The histograms are those of the 
simulation by model A and A-excitation. The dashed straight lines correspond to index of 
(a) 1.93 and (b) 2.66. Note that the muon antineutrino spectra are the same. 
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T 




[GeV] 

Fig. 14. — The 7r° gamma-ray spectra in the Galactic ridge obtained with Galprop with the 
parameter set, galdef 500180: the Galprop built-in p — p interaction model (dashed line); 
the present parameterized p — p interaction model (dash-dot line); the present parameterized 
model normalized to the dashed line in the energy band 0.7-0.8 GeV (solid line). 



- 29 - 



6. Application to Galactic Diffuse Gamma-Ray Emission 



We have replaced the 7r° production subroutine of Galprop (IStrong &: Moskalenkdll997l . 



2OOII ) with the present parameterized model and compared the Galactic diffuse gamma-ray 
spectra of tt" origin with that by the b uilt-in subroutine. A common parameter set, galdef 
500180 described in IStrong et al. J2004h . has been used in the two calculations. 



As shown in Figure [Ml the present parameterized model gives a flatter and smoother 
spectral energy distribution between = 0.3 — 2 GeV, smaller gamma-ray yield between 
0.5 — 1.3 G eV, and a hig h er po wer-law index between 1 — 5 GeV than the Galprop built-in 
model. In iKamae et al.l (120051 ). gamma-ray spectrum of the Galprop built-in model was 
compared with that of model A, after being normalized in the energy range E < 300 MeV. 
This normalization enhanced gamma-ray yield in the GeV range relative to sub-GeV range. 
The Galprop built-in p — p int eraction mode l has b een tuned to reproduce accelerator experi- 
ments better than model A of iKamae et al.l (120051 ) near the threshold. We have included the 
resonance contributions in the present model to improve this shortcoming near the thresh- 
old. Hence we have change the normalization point to the peak region in E'^dfinx/dE, 
E^ = 0.7 — 0.8 GeV, where gamma-rays from vr*^ decays are expected to dominate (see Fig- 
ure [HI dashed curve). The present model thus normalized gives ~ 20% higher gamma-ray 
yield than Galprop at 2 GeV: this is sub stantially lower than the difference (about 50%) 
shown in Figure 7 of iKamae et al.l (120051 ). The gamma-ray power-law index of the present 
mod el is harder by abou t 0.05 than that by the Galprop built-in model, just as the model 
A of iKamae et al.l (120051 ) predicts. 



7. Conclusion and Future Prospects 



We have presented the inclusive cross sections of stable secondary particles (7, e^, z/g, 
z/e, i/^, and Ufj^) produced by the p — p interaction in parameterized formulae. They facilitates 
computation of secondary particle spectra for arbitrary proton spectra as shown for Galactic 
diffuse gamma-ray spectrum (Figure [H]). Various effects that these secondary particles may 
have in astronomical environments can also be calculated at a higher precision. The formulae 
incorporate all important known features of the p — p interaction up to about Tp = 500 TeV 
and hence will also be useful in calculating background when searching for new phenomena. 

The parameterized model predicts all secondary particle spectra to have harder power- 
law indices than that of the incident proton and their inclusive cross sections to be larger 
than those expected fr om the old p — p interaction mode ls. When used to replace the p — p 



subroutine in Galprop (IStrong &: Moskalenkdll997l . l200ll ). the model gives a flatter spectral 
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energy density distribution between 0.3 — 5 GeV. The absolute gamma-ray yield predicted by 
the model is smaller than that by the Galprop model for < 1.5 GeV but higher for > 
1.5 GeV. If normalized near the peak in E'^dfiux/dE {E^ = 0.7 — 0.8 GeV), the parameterized 
model gives ~ 20% higher gamma-ray yield at 2 GeV than the model of Galpr op: our model 
with this normalization can account for ~ 20%, not ~ 50% as was claimed in Kamae et al. 



( I2OO5I ) , of the discrepancy between the diffuse Galactic Ridge gamma-ray spectrum observed 
by EGRET and r aodel predictions fo r the proton spectrum near the solar system (the power- 
law index ~ 2.7) (IHunter et al.lll997l ). We note the discrepancy w ill be reduced by inc lusion 
of the inverse Compton component (see, for example, a model by IStrong et al.l (120041 )). 



The present model also predicts more than e~ at energies higher wherever those 
produced hj p — p interaction become comparable in flux to primary e~ and e"*". The 
formulae and parameters given in the appendices are available as supplementary on-line 
material both in C language format and as C subroutine. We are currently parametrizing 
the angular distribution of gamma-rays relative to the incident proton direction. The results 
will be published elsewhere. 
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Table 1: Constants for equations [11 [21 [3] and [H 



ao 


= 0.1176 


bo = 11.34 


Co 


= 28.5 


do 


= 0.3522 


eo = 


5.922 


/o 


= 0.0834 


90 


= 0.0004257 




= 0.3829 


bi = 23.72 


Cl 


= -6.133 


do 


= 0.1530 


ei = 


1.632 


A 


= 9.5 


9i 


= 4.5 


02 


= 23.10 




C2 


= 1.464 


d2 


= 1.498 








= -5.5 


92 


= -7.0 


as 


= 6.454 








ds 


= 2.0 






h 


= 1.68 


93 


= 2.1 




= -5.764 








di 


= 30.0 






h 


= 3134 


94: 


= 503.5 




= -23.63 








^5 


= 3.155 














oe 


= 94.75 










= 1.042 















a7 = 0.02667 



Table 2: Kinematic limit parameters for the non-diffractive process. Proton kinetic energy 
Tp is in GeV. 



Particle 








7 


0.961ogio(Tp) 


15 


44 


e" 


0.961ogio(Tp) 


20 


45 


e+ 


0.941ogio(Tp) 


15 


47 




0.981ogio(Tp) 


15 


42 




0.981ogio(Tp) 


15 


40 




0.941ogio(Tp) 


20 


45 




0.981ogio(Tp) 


15 


40 
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Table 3: Parameters describing gamma-ray spectra for arbitrary proton energy. 

Parameter Foriiriilae as frmctions of the proton kinetic energy {y = log2[)(Tp)) in TeV. 
Non-diff . . equation EI 

-0.51187(1/ + 3.3) + 7.6179(1/ + 3.3)^ - 2.1332(y + 3.3)^ + 0.22184(j/ + 3.3)* 

-1.2592 ■ 10"^ + 1.4439 ■ lO'^exp {-0.29360(i/ + 3.4)) + 5.9363 ■ W-^l(y + 4.1485) + 2.2640 ■ lO""!/ - 3.3723 • 10" V 
-174.83 + 152.78 logi„ (1.5682(i/ + 3.4)) - 808.74/(i/ + 4.6157) 
0.81177 + 0.563851/ + 0.0040031i/^ - 0.0057658i/' + 0.000120572/* 

0.68631(1/ + 3.32) + 10.145(i/ + 3.32)^ - 4.6176(i/ + 3.32)^ + 0.86824(i/ + 3.32)* - 0.053741(i/ + 3.32)= 
9.0466 ■ 10"^ + 1.4539 ■ lO-'logn, (0.015204(i/ + 3.4)) + 1.3253 • 10-*/(i/ + 4.7171)^ - 4.1228 ■ lOSj + 2.2036 ■ \OSf 
-339.45 + 618.73 logi„ (0.31595(i/ + 3.9)) + 250.20/(t/ + 4.4395)^ 
-35.105 + 36.1671/ - 9.3575i/^ + 0.33717i/' 
0.17554 + 0.373001/ - 0.014938;/^ + 0.0032314i/3 + 0.0025579i/* 



r(y) 3.05 exp(-107((i/ + 3.25)7(1 + 8.08(1/ + 3.2.5)))^) for < 1.95 GeV 

1.01 for Tp > 1.95 GeV 

Diffraction, equation El 

To 60.142 tanh (-0.37555(t/ + 2.2)) - 5.9564(i/ + 0.-59913)^ + 6.0162 ■ IQ-'iy + 9.4773)* 

61 35.322 + 3.8026 tanh(-2.5979(i/ + 1.9)) - 2.1870 ■ 10-*(i/ + 369.13)^ 

62 -15.732 - 0.082064 tanh (-1.9621(i/ + 2.1)) + 2.3355 ■ 10-*(i/ + 252.43)^ 

63 -0.086827 + 0.37646 exp (-0.53053((i/ + 1.0444)/(1.0 + 0.27437(i/ + 1.0444)))^) 
Important: 60, .... 63 = for Tp < 5.52 GeV 

64 2.5982 + 0.39131(1/ + 2.95)^ - 0.0049693(i/ + 2.95)* + 0.94131 exp (-24.347(i/ + 2.45 - 0.19717(i/ + 2.45)^)^) 

65 0.11198 - 0.645821/ + OMlUy^ + 2.2853 exp (-0.0032432((i/ - 0.83562)/(1.0 + 0.33933(i/ - 0.83562)))^) 
bf, 1.7843 + 0.919141/ + 0.050118i/2 + 0.038096i/' - 0.027334i/* - 0.0035556i/= + 0.0025742i/« 
67 -0.19870 - 0.0710031/ + 0.019328i/^ - 0.28321 exp (-6.0516(i/ + 1.8441)^) 



A(1232), equationHa 

2.4316exp (-69.484((i/ + 3.1301)/(1.0 + 1.24921(y + 3.1301)))^) - 6.3003 - 9.5349/:y + 0.381211/^ 
56.872 + 40.6271/ + 7.75281/^ 

-5.4918 - 6.7872 tanh (4.7128(i/ + 2.1)) + 0.68048y 
-0.36414 + 0.0397771/ 

-0.72807 - 0.488281/ - 0.092876i/^ 

ResieOO), equationHa 



do 3.2433 exp (-57.133((i/ + 2.9507)/(1.0 + 1.2912(i/ + 2.9507)))^) - 1.0640 - 0.43925i/ 

<ii 16.901 + 5.95391/ - 2.12571/^ - 0.92057i/' 

<i2 -6.6638 - 7.5010 tanh (30.322(y + 2.1)) + 0.54662i/ 

rf;, -1.50648 - 0.872111/ - 0.17097i/^ 

rfj 0.42795 + 0.551361/ + 0.20707i/^ + 0.027552i/-' 



ao 

Qi 

02 
as 
ai 
as 
ao 
ai 
as 



Co 

Cl 
C2 

C3 
C4 
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Table 4: Parameter describing electron spectra for arbitrary proton ener^ 




Parameters 


Formulae as functions of the proton kinetic energy [y = logjQ(Tp)) in TeV. 




Non-diff. , equation [HI 


an 


-0.018639(1/ + 3.3) + 2.4315{y + 3.3)^ - 0.57719(!/ + 3.3)''^ + 0.063435(y + 3.3)-' 




ai 


7.1827 • 10-*' - 3.5067 ■ W-'^y + 1.3264 • lO^V - 3.3481 ■ lO^V + 2.3551 ■ 10-**:y^ + 3.4297 • 10"** 


y' 


02 


563.91 - 362.18 login (2-7187(?/ + 3.4)) - 2.8924 ■ lO'VCy + 7.9031)^ 




"3 


0.52684 + 0.57717y + 0.0045336j/2 - 0.0089066y^ 




Hi 


0.36108(y + 3.32) + 1.6963(y + 3.32)^ - 0.074456(y + 3.32)-' - 0.071455(y + 3.32)'' + 0.010473(.y 4 


-3.32)= 


15 


9.7387 • 10-= + 7.8573 ■ 10--= logjo (0.0036055(y + 4.3)) + 0.00024660/(!/ + 4.9390) - 3.8097 • Ur'y 


,2 


Qg 


-273.00 - 106.22 logm (0.34100(y + 3.4)) + 89.037)/ - 12.546?/2 




ar 


432.53 - 883.99 logio (0.19737(!/ + 3.9)) - 4.1938 ■ 107(?/ + 8.5518)^ 




as 


-0.12756 + 0.43478y - 0.0027797y2 - 0.0083074y=' 






3.63exp(-106((y + 3.26)/(l + 9.21(y + 3.26)))") - 0.182y - 0.175y" for < 15.6 GeV 
1.01 for Tp > 15.6 GeV 




Diffraction, equation IHl 




0.20463 tanh (-6.2370(y + 2.2)) - 0.16362(y + 1.6878)" + 3.5183 • 10-''(y + 9.6400)" 




h 


1.6537 + 3.8530 exp (-3.2027((y + 2.0154)/(1.0 + 0.62779(y + 2.0154)))") 




b, 

bi 
bs 


-10.722 - 0.082672 tanh (-1.8879(y + 2.1)) + 1.4895 • lO^^fe + 256.63)" 
-0.023752 - 0.51734 exp (-3.3087((y+ 1.9877)7(1.0 + 0.40300(y+ 1.9877)))") 
Important: 5o, . . . , 63 = for Tp < 5.52 GeV 

0.94921 + 0.12280(y + 2.9)" - 7.1585 • 10-''(y + 2.9)'' + 0.52130 log^ (y + 2.9) 
-4.2295 - 1.0025 tanh (9.0733(y + 1.9)) - 0.11452 * (y - 62.382) 




b. 


1.4862 + 0.99544y - 0.042763y" - 0.0040065y3 + 0.0057987y'' 




bj 


6.2629 + 6.9517 tanh (-0.36480(y + 2.1)) - 0.026033 * (y - 2.8542) 




ResieOO), equation in 


da 


0.37790 exp (-56.826((y + 2.9537)/(1.0 + 1.5221(y + 2.9537)))") - 0.059458 + 0.0096583y" 




di 


-5.5135 - 3.3988y 




ck 


-7.1209 - 7. 1850 tanh (30.801 (y + 2.1)) + 0.35108y 




d:. 


-6.7841 - 4.8385y - 0.91523y" 




di 


-134.03 - 139.63y - 48.316y" - 5.5526y^ 
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Table 5: Parameters describing positron spectra for arbitrary proton energy. 

Parameters Formulae as fmictions of the proton kinetic energy [y = log^Q(Tp)} in TeV. 
Non-diff . , equation [01 

-0.79606(y + 3.3) + 7.7496(y + 3.3)^ - 3.9326(y + 3.3)' + 0.80202(1/ + 3.3)* - 0.054994(j/ + 3.3)^ 
6.7943 ■ lO-** - 3.5345 ■ lO-'i/ + 6.0927 ■ 10" V + 2.0219 ■ lO-'y' + 5.1005 ■ 10" V - 4.2622 ■ 10"")/= 
44.827 - 81.378 logio {0.027733(y + 3.5)) - 1.3886 ■ 10V(:y + 8.4417) 
0.52010 + 0.593361/ + 0.012032i/2 - 0.0064242i/' 

2.1361(!/ + 3.32) + 1.8514(i/ + 3.32)^ - 0.47872(j/ + 3.32)' + 0.0032043(y + 3.32)'' + 0.0082955(?/ + 3.32)^ 
1.0845 ■ 10-<' + 1.4336 ■ 10"" logn, {0.0077255(i/ + 4.3)) + 1.3018 • lO'Vly + 4.8188)^ + 9.3601 ■ 10"**)/ 
-267.74 + 14.175 logio (0.35391(i/ + 3.4)) + 64.669/(i/ - 7.7036)=^ 
138.26 - 539.84 logio {0.12467(i/ + 3.9)) - 1.9869 ■ 107(i/ + 7.6884)=^ + 1.0675t/^ 
-0.14707 + 0.401351/ + 0.0039899i/' - 0.0016602i/'' 



r(y) 2.22exp(-98.9((2/ + 3.25)/(l + 10.4(y + 3.25)))^) for < 5.52 GeV 

1.0 for Tp > 5.52 GeV 
Diffraction, equation El 

Ta 29.192 tanh (-0.37879(j/ + 2.2)) - 3.2196(i/ + 0.67500)=^ + 0.0036687(i/ + 9.0824)" 

61 -142.97 + 147.86exp (-0.37194((i/ + 1.8781)/(1.0 + 3.8389(j/ + 1.8781)))^) 

62 -14.487 - 4.2223 tanh (-13.546(j/ + 2.2)) + 1.6988 ■ IQ-'^iy + 234.65)^ 

63 -0.0036974 - 0.41976exp (-6.1527((i/ + 1.8194)/{1.0 + 0.99946(1/ + 1.8194)))^) 
Important: bg, . . . ,63 = loi < 11.05 GeV 

64 1.8108 + 0.18545(1/ + 2.9)^ - 0.0020049(1/ + 2.9)* + 0.85084exp (-14.987(.t + 2.29 - 0.18967(i + 2.29)^)^) 

65 2.0404 - 0.51548 tanh (2.2758{)/ + 1.9)) - 0.035009/(1/ - 6.6555) 
bf, 1.5258 + 1.01321/ - 0.0643881/^ - 0.0040209i/' - 0.0082772y* 

67 3.0551 + 3.5240 tanh (-0.36739(y + 2.1)) - 0.13382 * (y - 2.7718) 



A(1232), equation m 

2.9841 exp (-67.857(()/ + 3.1272)/(1.0 + 0.22831(j/ + 3.1272)))=^) - 6.5855 - 9.6984/y + 0.41256i/2 
6.8276 + 5.2236J/ + 1.4630t/^ 

-6.0291 - 6.4581 tanh(5.0830(y + 2.1)) + 0.46352i/ 
0.59300 + 0.360931/ 

0.77368 + 0.447761/ + 0.056409;/^ 

ResieOO), equation 113 



do 1.9186 exp (-56.544({i/ + 2.9485)/(1.0 + 1.2892(i/ + 2.9485)))^) - 0.23720 + 0.041315^^ 

<ii -4.9866 - 3.1435t/ 

<i2 -7.0550 - 7.2165 tanh(31.033(t/ + 2.1)) + 0.38541i/ 

rf;, -2.8915 - 2.14951/ - 0.45006i/^ 

rfj -1.2970 - 0.139471/ - 0.41197i/^ - 0.10641!/'' 



ao 
ai 
02 
as 
ai 
as 
ao 
ai 
as 



Co 

Cl 
C2 

C3 
C4 
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Table 6: Parameters describing electron neutrino spectra for arbitrary proton energy. 

Parameters Formulae as functions of the proton kinetic energy {y = logj,){Tp)) in TeV. 



Non-diff., equation 





0.0074087 + 2.9161(7/ + 3.31) + 0.99061(y + 3.31)^ - 0.28694(j/ + 3.31)=' + 0.038799(y + 3.31)" 




Ql 


-3.2480 • 10-» + 7.1944 ■ lO^^^ exp (-0.21814(j/ + 3.4)) + 2.0467 ■ Kr^iy + 4.1640) + 5.6954 ■ 


lO-'V - 3.4105 ■ 10-^y^ 


12 


-230.50 + 58.802y - 9.9393y2 + 1.2473^^ - 0.26322y'' 




a-i 


0.45064 + 0.56930y + 0.012428j/2 - 0.0070889y=* 




a,! 


-0.011883 + 1.7992(j/ + 3.32) + 3.b2U{y + 3.32)^ - 1.7478(y + 3.32)=* + 0.32077(|/ + 3.32) ' - 


0.017667(y + 3.32)=^ 


Qs 


-1.6238 • 10-' + 1.8116 ■ 10-^ exp (-0.30111(j/ + 3.4)) + 9.6112 ■ W-y(y + 4.8229)^ 




as 


-261.30 - 43.351 log^ (0.35298(y + 3.4)) + 70.925/(y - 8.7147)=^ 




a? 


184.45 - 1473.6/(y + 6.8788) - 4.0536^=^ 




as 


-0.24019 + 0.38504y + 0.00968691/^ - 0.0015046^' 





r{y) 0.329 exp (-247((y + 3.26)/(l + 6.56(y + 3.26)))") - 0.957y - 0.229?/=' for T, < 7.81 GeV 



1.0 for > 7.81 GeV 

Diffraction, equation 01 

To 53.809 tanlr {-0.41421(y + 2.2)) - 6.7538(y + 0.76010)=^ + 0.0088080(j/ + 8.5075)" 

61 -50.211 + 55.131 exp (1.3651((y + 1.8901)/(1.0 + 4.4440(y + 1.8901)))") 

62 -17.231 + 0.041100 tanh {7.9638(y + 1.9)) - 0.055449y + 2.5866 ■ 10-''(y + 250.68)" 

63 12.335 - 12.893 exp (-1.4412((y + 1.8998)/(1.0 + 5.5969(y + 1.8998)))") 
Important: 5o, . . . , 63 = for Tp < 11.05 GeV 

64 1.3-558 + 0.46601(y + 2.95) + 0.052978(y + 2.2)" + 0.79575 exp (-5.4007(y + 2.2 + 4.6121(2: + 2.2)")") 

65 1.8756 - 0.42169 tanli (1.6100(y + 1.9)) - 0.051026 * (y - 3.9573) 
bn 1.5016 + 1.0118y - 0.072787y" - 0.0038858y=' + 0.0093650y" 

67 4.9735 + 5.5674 tanh (-0.36249(?7 + 2.1)) - 0.20660 * (y - 2.8604) 

Ari2321. eouationini 

~o 2.8290 exp (-71.339((y + 3.1282)/(1.0 + 0.48420(y + 3.1282)))") - 9.6339 - 15.733/y + 0.52413y" 

Cl -24.571 - 15.831y- 2. 1200y" 

C2 -5.9593 - 6.4695 tanh (4.7225(y + 2.1)) + 0.50003y 

C3 0.26022 + 0.24545y 

C4 0.076498 + 0.061678y + 0.0040028y" 

Resl600), eriuationlH 

dji 1.7951 exp (-57.260((y + 2.9509)/(1.0 + 1.4101(y + 2.9509)))") - 0.58604 - 0.23868y 

di -2.6395 - 1.5105y + 0.22174y" 

d2 -7.0512 - 7.1970tanh(31.074(y + 2.1)) + 0.39007y 

d-i -1.4271 - 1.0399y-0.24179y" 

di 0.74875 + 0.63616y + 0.17396y" + 0.017636y=' 
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Table 7: Parameters describing electron anti-neutrino spectra for arbitrary proton energy. 

Parameters Formulae as functions of the proton kinetic energy {y = logio{I]i)) in TeV. 



Non-diff. , equation [U] 





O.OOf 3f f3 + 0.36538(y + 3.3f ) + 1.5178(y + 3.31)^ - 0.20668(|/ + 3.3f + 0.024255(y + 3.3f )' 




ai 


-4.7833 • f 0-'' + 4.5837 ■ IQ-'" exp (-0.42980(i/ + 3.4)) + 6.1559 ■ 10-y{y + 4.1731) + 1.1928 ■ 10" 




0.2 


-245.22 + 73.223)/ - 19.652j/2 + 0.083138i/ + 0.71561y'' 




a-i 


0.45232 + 0.52934y + 0.010078y" - 0.0017092^^* 




ai 


-0.0025734 + 0.38424(y + 3.32) + 1.5517(y + 3.32)^ + 0.17336()/ + 3.32)'' - 0.17160(y + 3.32)'* + 


0.021059(y + 3.32)» 


15 


4.7673 • 10-"* + 5.4936 ■ IQ-'* log (0.0067905(y + 4.3)) + 0.00020740/(j/ + 4.9772) 




ae 


-270.30 - 114.471ogio (0.34352(y + 3.4)) + 80.085)/ - 7.9240j/" 




ai 


3272.9 - 2.9161 ■ 10V(y + 87.847) - 6.2330y" 




a.% 


-0.17787 + 0.36771y - 0.025397y^ + 0.0019238y=' + 0.0032725^" 




r{y) 


2.67 exp (-45.7((y + 3.27)/(l + 6.59{y + 3.27)))^) - 0.301)/ - 0.208y" for Tp < 15.6 GeV 
1.0 for Tp > 15.6 GeV 




Diffraction, equation 01 


60 


41.307 tanli (-0.37411(y + 2.2)) - 4.1223(y + 0.55505)" + 0.0042652(y + 9.2685)" 




61 


-132.50 + 142.12 exp (-8.0289((y+ 1.9196)7(1.0+ 11. 530(y+ 1.9196)))") 




b2 


-17.223 + 0.011285 tanh(-69.746(y + 1.9)) - 0.048233y + 2.5881 • 10-''(y + 250.77)" 




h 


8.1991 - 9.6437 exp (-45.261((y + 1.9292)/(1.0 + 16.682(y + 1.9292)))") 
Important: 60, ■ ■ ■ , 63 = for Tp < 11.05 GeV 




h 


0.55919 + 0.36647{y + 2.95)" + 0.056194(y + 2.95)* + 0.49957exp (-5.5317(y + 2.2 + 0.43867(y + 2.2)")") 


h 


1.2544 - 0.52362 tanh (2.7638(y + 1.9)) + 0.055837 * (y - 17.638) 






1.4788 + 1.0278y - 0.092852y" - 0.0062734y3 + 0.011920y'' 




67 


5.1651 + 5.7398 tanh (-0.37356(y + 2.1)) - 0.22234 * (y - 2.7889) 




ResieOO), 


equation [121 




do 


0.36459 exp (-58.210((y + 2.9537)/(1.0 + 1.4320(y + 2.9537)))") - 0.11283 - 0.046244y 




rfi 


-9.5066 - 5.4655y - 0.31769y" 




d2 


-7.1831 - 7.1551ian/i(30.354(y + 2.1)) + 0.33757y 




d:. 


2.7938 + 1.6992y + 0.20161y 




di 


0.61878 + 0.62371y + 0.18913y" + 0.019118y'' 
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Table 8: Parameters describing muon neutrino spectra for arbitrary proton energy. 

Parameters Formulae as functions of the proton kinetic energy {y = logj(){Tp)) in TeV. 



Non-diff. , equation 





-0.63611()/ + 3.3) + 9.9015(y + 3.3)^ - 4.5897(y + 3.3)-' + 0.91778{y + 3.3)* - 0.060724(y + 3.3)-* 


ai 


6.8700 ■ 10-'' - 2.8245 ■ lO"''?; + 7.6032 ■ lO^V - 3.2953 ■ lO^V + 7.4292 ■ ur'^y'^ 


0.2 


-240.46 + 58.405y - 9.8556y^ + 3.1401^^ - 0.88932y^ 


a-i 


0.49935 + 0.60919y + 0.0024963!/2 - 0.0099910^^ 


ai 


2.5094(j/ + 3.32) + 4.1350(|/ + 3.32)^ - 0.89534()/ + 3.32)=' - 0.0027577(y + 3.32)'' + 0.014511{y + 3.32)^ 


as 


8.2046 ■ 10-^ + 1.4085 ■ lO"" logi„ (0.016793(!/ + 4.3)) + 0.00013340/(y + 4.7136)^ 


as 


-267.55 - 0.21018 logio (0.35217(y + 3.9)) + 69.586?/ - 9.9930^'' 


a? 


2742.8 + 222.01 logi„ (9.7401(y + 3.9)) - 4772.5/(y + 19.773) - 6.1001y" 


as 


-0.11857 + 0.390721/ - 0.037813i/ + 0.0022265y=' + 0.0046931)/'' 


r{y) 


2.23 exp (-93.4((|/ + 3.25)/(l + 8.38()/ + 3.25)))^^) - 0.376|/ - 0.1211/^ for < 15.6 GeV 




1.0 for Tp > 15.6 GeV 


Diffraction, equation IHl 


60 


64.682 tanh (-0.34313(i/ + 2.2)) - 5.5955(;/ + 0.44754)^ + 0.0050117(?/ + 9.9165)" 




-7.6016 + 3.0427 • 10" exp (-1.0134 ■ 10"((y + 2.3066)/(1.0 + 41.612(y + 2.3066)))") 


62 


-1.4978 - 0.58163 tanli(-0.36488(j/ + 1.9)) + 0.031825(j/ + 2.8097) + 0.022796(y - 1.8861)" 


63 


-0.0061483 - 65.799 exp (-4.8239((y + 3.8835)/(1.0 + 0.53343()/ + 3.8835)))") 




Important: feo, . . . , 63 = for Tp < 11.05 GeV 


h 


2.8009 + 0.35341(y + 2.95)" - 0.0039779(y + 2.95)" + 1.3012 exp (-10.592(y + 2.2 - 0.19149(y + 2.2)")") 


h 


1.8016 - 0.69847 tanli(2.8627(y + 1.9)) - 0.015722 * (y - 45.4) 


K 


1.4617 + 1.0167y - 0.078617y" - 0.0038336y'' + 0.010141y" 




3.5599 + 4.0041 tanh (-0.41889(j/ + 2.1)) - 0.18182 * {y - 2.4209) 


A(1232), equation ini 


Co 


3.6052 exp (-60.914((y + 3.1278)/(1.0 - 0.19497(y + 3.1278)))") - 0.92514 + 2.1315/y + 0.23548y" 


Cl 


95.310 + 70.497y+ 13.636y" 


C2 


-6.2158 - 6.2939 tanh(21.592(j/ + 2.1)) + 0.37440y 


C3 


2.7485 + 1.1692y 


C4 


-2.7568 - 1.8461y - 0.31376y" 


Resf600), eriuation[T3 


do 


2.5489 exp (-58.488((y + 2.9509)/(1.0 + 1.3154(y + 2.9509)))") - 0.83039 - 0.34412y 


d. 


88.173 + 65.148y + 12.585y" 


d2 


-7.0962 - 7.1690 tanh(30.890(y + 2.1)) + 0.38032y 


d:. 


-4.1440 - 3.2717y - 0.70537y" 


di 


2.2624 + 1.1806y - 0.0043450y" - 0.043020r' 
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Table 9: Parameters describing muon anti-neutrino spectra for arbitrary proton energy. 

Parameters Formulae as functions of the proton kinetie energy {y = logj(,(Tp)) in TeV. 



Non-diff . , equation [B] 



fflO 


-1.5243(2/ + 3.3) + 10.107(2/ + 3.3)^ - 4.3126(?/ + 3.3)^ + 0.80081(;y + 3.3)* - 0.048724(2/ + 3.3)= 


ai 


-2.6297 • 10-5 + 9.3858 ■ 10"= exp (-0.32384(2/ + 3.4)) + 7.7821 • 10-V(;y + 4.0560) + 7.6149 • 10-'^;y - 8.4091 ■ 10" V 




-223.62 + 59.3742/ - 5.73562/^ + 1.98152/^ - lM78y* 


ffls 


0.50807 + 0.60221?/ + 0.00341202/^ - 0.0111392/=* 




2.6483(2/ + 3.32) + 4.4585(2/ + 3.32)^ - 1.2744(2/ + 3.32)=' - 0.11659(2/ + 3.32)'' + 0.0030477(y + 3.32)= 


as 


9.1101 ■ 10-' + 1.3880 • 10-'' logic (0.016998(2/ + 4.3))1.2744 • W-*/{y + 4.7707)^ 


ae 


-272.11 - 53.477 logio (0.35531(t/ + 3.9)) + 56.041/(2/ - 6.0876)^ 


a? 


6431.8 + 893.92 logy, (5.713 ■ 10'^[y + 3.9)) + 2 103.6/(2/ + 5.67 40) - 6.11252/^ 


08 


-0.11120 + 0.381442/ - 0.0401282/^ + 0.00474842/=^ + 0.00547072/'' 



r(2/) 2.56exp (-107((2/ + 3.25)/(l + 8.34(2/ + 3.25)))2) - 0.3852/ - 0.1252/^ for Tp < 15.6 GeV 
1.0 for Tp > 15.6 GeV 



Diffraction, equation El 

bii 70.430 tanh (-0.35816(2/ + 2.2)) - 6.6796(2/ + 0.52273)^ + 0.0065659(2; + 9.5266)'' 

61 -8.1145 + 7686.0exp (4.4046 • 10\{y + 2.2190)/(1.0 + 81.105(2/ + 2.2190)))^) 

62 -1.3095 + 0.071270 tanh (-0.0075463(2/ + 1.9)) + 0.067759(2/ + 5.3433) - 0.0044205('(/ - 1.8683)= 

63 0.082149 - 2190.1 exp (-533.75((2/ + 2.8363)/(1.0 + 7.0976(2/ + 2.8363)))=) 
Important: 6,), . . . , 63 = for Tp < 11.05 GeV 

64 2.7540 + 0.33859(2/ + 2.95)= - 0.0034274(2/ + 2.95)* + 1.1679 exp (-10.408(2/ + 2.2 - 0.18922(2/ + 2.2)=)=) 

65 2.1817 - 0.59584 tanh (2.7054(2/ + 1.9)) - 0.010909 * (y - 14.9) 
be 1.4591 + 1.02752/ - 0.0749492/= - 0.00603962/=* + 0.00 9 75682/" 

67 3.7609 + 4.2843tanh (-0.37148(2/ + 2.1)) - 0.16479* (2/ - 2.7653) 



A(1232), equationlH 



Co 


2.8262 exp (-62.894((2/ + 3.1250)/(1.0 - 


0.47567(2/ + 3.1250)))=) + 5.6845 + 13.409/2/ - 0.0972962/= 


Cl 


16.721 + 11.750:y + 2.4637;;/= 




C2 


-6.0557 - 6.3378 tanh(21.984(2/ + 2.1)) 


+ 0.431732/ 


C3 


0.37009 + 0.277062/ 




C4 


0.047507 + 0.0615702/ + 0.0070117y= 





ReslOOO), equation O 



do 2.2400 exp (-57.159((y + 2.9492)/(1.0 + 1.2994(2/ + 2.9492)))=) - 0.66521 - 0.27554y 

di -7.0650 - 4.27732/ - 0.176482/= 

di -7.0410 - 7.1977 tanh(31.095(2/ + 2.1)) + 0.40238;y 

dg -1.2354 - 0.875812/ - 0.20829= 

d4 -0.11395 + 0.344182/ + 0.27103;y= + 0.050248:y=* 



